A cost-effective blood DNA methylation-based age estimation method in domestic cats, Tsushima leopard cats (Prionailurus bengalensis euptilurus) and Panthera species, using targeted bisulphite sequencing and machine learning models.

Individual age can be used to design more efficient and suitable management plans in both in situ and ex situ conservation programmes for targeted wildlife species. DNA methylation is a promising marker of epigenetic ageing that can accurately estimate age from small amounts of biological material, which can be collected in a minimally invasive manner. In this study, we sequenced five targeted genetic regions and used 8-23 selected CpG sites to build age estimation models using machine learning methods at only about $3-7 per sample. Blood samples of seven Felidae species were used, ranging from small to big, and domestic to endangered species: domestic cats (Felis catus, 139 samples), Tsushima leopard cats (Prionailurus bengalensis euptilurus, 84 samples) and five Panthera species (96 samples). The models achieved satisfactory accuracy, with the mean absolute error of the most accurate models recorded at 1.966, 1.348 and 1.552 years in domestic cats, Tsushima leopard cats and Panthera spp. respectively. We developed the models in domestic cats and Tsushima leopard cats, which were applicable to individuals regardless of health conditions; therefore, these models are applicable to samples collected from individuals with diverse characteristics, which is often the case in conservation. We also showed the possibility of developing universal age estimation models for the five Panthera spp. using only two of the five genetic regions. We do not recommend building a common age estimation model for all the target species using our markers, because of the degraded performance of models that included all species.

Analysis of Crop Consumption Using Scatological Samples from the Red-Crowned Crane Grus japonensis in Eastern Hokkaido, Japan.

Total DNA extracts from the intestinal contents of 60 flying red-crowned cranes (juveniles, subadults and adults) found dead in 2006-2021, and the feces of 25 chicks collected in June and July of 2016-2018, were used for PCR reactions with primers specific for 16 crops, followed by high-throughput sequencing. The most predominant crop detected was corn in adult and subadult cranes (61.7%). Other grains (barley, wheat, soybean) (5.0-8.3%) and vegetables (tomatoes, Chinese cabbage, etc.) (1.7-6.7%) were also detected in flying cranes. Surprisingly, some of the detected crops were not grown in the Kushiro and Nemuro regions. There was no significant difference in crop intake status in winter and that in other seasons for most of the crops. Corn (28.0%), soybeans (8.0%), wheat and beet (4.0%) were detected in crane chicks in summer, though the detection rates were generally lower than those in flying cranes. Alfalfa, which is not grown in eastern Hokkaido but is used in some cattle feed, was detected in some cranes. Rice, buckwheat, adzuki beans, common beans, potatoes and carrots were not detected at any life stage, indicating the preferences of red-crowned cranes. The results suggest that red-crowned cranes in Hokkaido are dependent on dairy farmers for their feed supply.

Cultured fibroblasts of the Okinawa rail present delayed innate immune response compared to that of chicken.

The Okinawa rail is endemic to Okinawa Island and is categorized as an endangered animal. In this study, we focused on innate immunity because it is the first line of host defense. In particular, signals recognizing foreign RNA (e.g., viruses) are important for host defense because they activate the host immune system. The retinoic acid-inducible gene I (RIG-I)-like receptor (RLR) families (RIG-I, MDA5, and LGP2) are sensors that activate innate immunity. Therefore, we analyzed these functions in the Okinawa rail using genomic and cellular analyses of fibroblasts. Fibroblasts can be obtained from dead individuals, allowing these cells to be obtained from dead individuals, which is particularly useful for endangered species. The MDA5 gene of Okinawa rail was sequenced using the Sanger method following PCR amplification and extraction of the amplified sequence from agarose gel. Additionally, mRNA expression analysis of cultured fibroblasts exposed to poly I:C was done. The MDA5 gene was found to be a mutated nonfunctional gene in the Okinawa rail. The mRNA expression rates of inflammatory cytokine genes type I IFN, and Mx1 were slower in Okinawa rail than in chicken cultured fibroblasts. Similar to the mRNA expression results, cell number and live cell ratio also slowly decreased in the Okinawa rail compared with chicken cultured fibroblasts, indicating that the innate immune reaction differs between chicken and the Okinawa rail. To the best of our knowledge, this is the first experimental evaluation of the loss of function of the Okinawa rail innate immune genes. In conclusion, our results provide a basis for conservation strategies for the endangered Okinawa rail.

The complete mitochondrial genome of water flea Ceriodaphnia dubia (Crustacea: Cladocera) NIES strain.

Water flea Ceriodaphnia dubia has been widely used for risk assessments of chemicals and environmental contamination. In this study, the complete mitochondrial genome (mitogenome) of this species NIES strain was determined using short-read high throughput and long-read sequencing technologies. The mitogenome of C. dubia was 15,170 bp in length and consisted of 13 protein-coding genes (PCGs), 2 ribosomal RNAs (rRNAs), and 22 transfer RNAs (tRNAs). The gene order was identical to the pattern conserved across crustaceans. The complete mitogenome of the NIES strain will serve as genetical reference in ecological risk assessments in Japan, as well as resources for future phylogenetical studies using cladocerans.

Draft genome sequence of Pelosinus sp. IPA-1, a bacterium isolated from arsenic-contaminated soil in Japan.

Pelosinus sp. strain IPA-1 is a bacterium isolated from arsenic-contaminated soil in Japan. We here report the draft genome sequence of strain IPA-1.

Draft genome sequences of Pantoea sp. strains QMID1-QMID4 isolated from the midgut of Japanese honey bee (Apis cerana japonica).

We report the draft genome sequences of Pantoea sp. strains QMID1-QMID4 that were recovered from the midgut of Japanese honey bee (Apis cerana japonica). The strains possess the carotenoid biosynthetic gene cluster. The genome information expands our knowledge of their potential use as probiotics and/or prebiotics in honey bees.

First report of the mitogenome of the invasive reef-building polychaete Ficopomatus enigmaticus (Annelida: Serpulidae) and a cryptic lineage from the Japanese Archipelago.

BACKGROUND: The mitochondrial genomes (mitogenomes) of the family Serpulidae are characterized by a high nucleotide sequence divergence and a significant number of gene order rearrangements compared with other families within the phylum Annelida. However, only two of 50 genera of serpulids have mitogenomes already sequenced. In this study, we report the first sequencing and assembly of the complete mitogenome of Ficopomatus, thus providing further knowledge on mitochondrial gene sequences of Serpulidae. METHODS AND RESULTS: A mitogenome of the invasive reef-building polychaete Ficopomatus enigmaticus was amplified by long PCR and sequenced using the Illumina MiSeq System. It comprised 15,853 bp and consisted of 12 protein-coding genes (atp8 was not found), 23 tRNA, and two rRNA genes. The AT and GC skew values were infrequent when compared to annelid mitogenomes but similar to other serpulids sequenced to date (i.e., Spirobranchus and Hydroides). The mitochondrial gene order of F. enigmaticus was highly rearranged compared to other serpulids. To amplify 16S rRNA gene sequences, we developed a 16S rRNA primer set by modifying the universal primer set 16SarL/16SbrH. We detected the 16S rRNA sequence of F. enigmaticus deposited in GenBank erroneously characterized as of serpulid origin. We reported for the first time the presence of two lineages of F. enigmaticus in Japan, which have already been identified in California, Australia, and the Mediterranean. CONCLUSIONS: The first mitochondrial genome of F. enigmaticus showed a unique gene order rearrangement, corroborating the remarkable diversity in the previously reported mitogenomes of other serpulid species. The presence of the two lineages of F. enigmaticus identified for the first time in Japan represents another case of cryptic invasion. The first 16S rRNA gene sequences of F. enigmaticus obtained in the present study can be used as reference sequences in future DNA metabarcoding studies.

Organoids containing neural-like cells derived from chicken iPSCs respond to poly:IC through the RLR family.

There is still much room for development in pluripotent stem cell research on avian species compared to human stem cell studies. Neural cells are useful for the evaluation of risk assessment of infectious diseases since several avian species die of encephalitis derived from infectious diseases. In this study, we attempted to develop induced pluripotent stem cells (iPSCs) technology for avian species by forming organoids containing neural-like cells. In our previous study, we established two types iPSCs from chicken somatic cells, the first is iPSCs with PB-R6F reprogramming vector and the second is iPSCs with PB-TAD-7F reprogramming vector. In this study, we first compared the nature of these two cell types using RNA-seq analysis. The total gene expression of iPSCs with PB-TAD-7F was closer to that of chicken ESCs than that of iPSCs with PB-R6F; therefore, we used iPSCs with PB-TAD-7F to form organoids containing neural-like cells. We successfully established organoids containing neural-like cells from iPSCs using PB-TAD-7F. Furthermore, our organoids responded to poly:IC through the RIG-I-like receptor (RLR) family. In this study, we developed iPSCs technology for avian species via organoid formation. In the future, organoids containing neural-like cells from avian iPSCs can develop as a new evaluation tool for infectious disease risk in avian species, including endangered avian species.

Molecular Phylogeny of Thoracotreme Crabs Including Nine Newly Determined Mitochondrial Genomes.

Mitochondrial genomes are used widely for the molecular phylogenetic analysis of animals. Although phylogenetic analyses based on the mitogenomes of brachyurans often yield well-resolved phylogenies, most interfamilial phylogenetic relationships in Thoracotremata remain unclear. We determined nine new mitogenomes of Thoracotremata, including mitogenomes of Camptandriidae (Deiratonotus japonicus), Dotillidae (Ilyoplax integra, Ilyoplax pusilla, and Tmethypocoelis choreutes), Macrophthalmidae (Ilyograpsus nodulosus), Pinnotheridae (Arcotheres sp. and Indopinnixa haematosticta), Plagusiidae (Guinusia dentipes), and Percnidae (Percnon planissimum). Interestingly, Percnon planissimum (Percnidae) was found to possess ≥ 19 repeated sequences in the control region. The gene orders of Il. nodulosus, Arcotheres sp., and In. haematosticta were revealed to be unique among thoracotreme crabs. Although the results of Bayesian and maximum likelihood (ML) phylogenetic analyses of three datasets were incongruent, highly supported clades (PP ≥ 0.99 or BS ≥ 99%) were not contradictory among the analyses. All analyses suggested the paraphyly of Grapsoidea and Ocypodoidea, corroborating the findings of previous studies based on molecular phylogenies of thoracotreme crabs. The phylogenetic positions of symbiotic thoracotreme crabs, Pinnotheridae and Cryptochiridae, were highly supported (Pinnotheridae + Ocypodidae and Cryptochiridae + Grapsidae, respectively) for the Bayesian analyses but not for the ML analyses. Analyses of more thoracotreme species' mitogenome sequences in additional studies will further strengthen the framework for thoracotreme evolution.

The dynamics of the microbiome in Ixodidae are shaped by tick ontogeny and pathogens in Sarawak, Malaysian Borneo.

Tick-borne diseases have recently been considered a potential emerging public health threat in Malaysia; however, fundamental studies into tick-borne pathogens and microbiome appear limited. In this study, six tick species (Ixodes granulatus, Haemaphysalis hystricis, Haemaphysalis shimoga, Dermacentor compactus, Dermacentor steini and Dermacentor atrosignatus) collected from two primary forests and an oil palm plantation in Sarawak, Malaysian Borneo, were used for microbiome analysis targeting bacterial 16S rDNA using next-generation sequencing (NGS). In addition, bacterial species were further characterized in conventional PCRs to identify potential pathogens. Sequences generated from NGS were first filtered with the Decontam package in R before subsequent microbial diversity analyses. Alpha and beta analyses revealed that the genus Dermacentor had the highest microbial diversity, and H. shimoga significantly differed in microbial composition from other tick species. Alpha and beta diversities were also significantly different between developmental stages of H. shimoga. Furthermore, we observed that some bacterial groups were significantly more abundant in certain tick species and developmental stages of H. shimoga. We tested the relative abundances using pairwise linear discriminant analysis effect size (LEfSe), which also revealed significant microbial composition differences between Borrelia-positive and Borrelia-negative I. granulatus ticks. Finally, pathogenic and potentially pathogenic bacteria circulating in different tick species, such as Rickettsia heilongjiangensis, Ehrlichia sp., Anaplasma sp. and Bartonella spp. were characterized by PCR and sequencing. Moreover, Coxiella and Francisella-like potential symbionts were identified from H. shimoga and D. steini, respectively. More studies are required to unravel the factors associated with the variations observed in this study.

Phytocyanin-encoding genes confer enhanced ozone tolerance in Arabidopsis thaliana.

Ozone is a phytotoxic air pollutant that has various damaging effects on plants, including chlorosis and growth inhibition. Although various physiological and genetic studies have elucidated some of the mechanisms underlying plant ozone sensitivity and lesion development, our understanding of plant response to this gas remains incomplete. Here, we show evidence for the involvement of certain apoplastic proteins called phytocyanins, such as AtUC5, that protect against ozone damage. Two representative ozone-inducible responses, chlorosis and stomatal closure, were suppressed in AtUC5-overexpressing plants. Analysis of transgenic plants expressing a chimeric protein composed of AtUC5 fused to green fluorescent protein indicated that this fusion protein localises to the apoplast of plant cells where it appears to suppress early responses to ozone damage such as generation or signalling of reactive oxygen species. Moreover, yeast two-hybrid analyses suggest that AtUC5 may physically interact with stress-related proteins such as copper amine oxidase and late embryogenesis abundant protein-like protein. In addition to AtUC5, other examined phytocyanins such as AtUC6 and AtSC3 could confer ozone tolerance to plants when overexpressed in A. thaliana, suggesting that these proteins act together to protect plants against oxidative stress factors.

Induced pluripotent stem cells of endangered avian species.

The number of endangered avian-related species increase in Japan recently. The application of new technologies, such as induced pluripotent stem cells (iPSCs), may contribute to the recovery of the decreasing numbers of endangered animals and conservation of genetic resources. We established novel iPSCs from three endangered avian species (Okinawa rail, Japanese ptarmigan, and Blakiston's fish owl) with seven reprogramming factors (M3O, Sox2, Klf4, c-Myc, Nanog, Lin28, and Klf2). The iPSCs are pluripotency markers and express pluripotency-related genes and differentiated into three germ layers in vivo and in vitro. These three endangered avian iPSCs displayed different cellular characteristics even though the same reprogramming factors use. Japanese ptarmigan-derived iPSCs have different biological characteristics from those observed in other avian-derived iPSCs. Japanese ptarmigan iPSCs contributed to chimeras part in chicken embryos. To the best of our knowledge, our findings provide the first evidence of the potential value of iPSCs as a resource for endangered avian species conservation.

Detection of a Babesia sp. genotype closely related to marsupial-associated Babesia spp. in male Haemaphysalis shimoga from Sarawak, Malaysian Borneo.

In this study, Babesia screening was conducted in 55 rodents and 160 tick samples collected from primary forests and an oil palm plantation in Sarawak, Malaysian Borneo. PCR targeting the 18S ribosomal DNA revealed the presence of Babesia spp. DNA detected in two questing male Haemaphysalis shimoga ticks collected from the oil palm plantation. Sequence analysis revealed that both sequences were identical and had 98.6% identity to a Babesia macropus sequence obtained from Eastern grey kangaroos (Macropus giganteus) in Australia. Phylogenetic tree revealed clustering with marsupial-associated Babesia spp. in the Babesia sensu stricto clade. Whether or not H. shimoga is the competent vector and the importance of the Babesia sp. detected in this study warrants more investigation.

Conventional Gel Electrophoresis-Resolvable Insertion/Deletion Markers for Individual Identification and Analysis of Population Genetics in Red-Crowned Cranes in Eastern Hokkaido, Japan.

Red-crowned crane Grus japonensis is an endangered species in two separate populations: the mainland population in the Eurasian continent and the island population in eastern Hokkaido, Japan. We found 11 insertion/deletion (InDel) markers in the genome of the red-crowned crane and designed primer sets across these InDels that can be analyzed with conventional agarose gel electrophoresis. Sixty-six samples of whole blood and skeletal muscle obtained from red-crowned cranes, including 12 families in eastern Hokkaido from 1994 to 2021, showed different patterns in gel images of 11 InDel PCR reactions except for two pairs. The combined non-exclusion probability of the 11 markers indicates that individuals can be determined with a probability of 99.9%. In 39 non-relative chicks, the expected heterozygosity (He) was 0.316, suggesting low genetic diversity. This might not be caused by high levels of inbreeding since the average FIS was not significantly different from zero (0.095, p = 0.075). The results suggest that the 11 InDel primer sets can be used for fairly accurate individual identification as well as genetic population analyses in red-crowned cranes in the island population.

Comparison of the Genetic Diversity of the Captive and Wild Populations of the Tsushima Leopard Cat Using a GRAS-Di Analysis.

The Tsushima leopard cat (Prionailurus bengalensis euptilurus) (TLC) is a regional population of the Amur leopard cat (P. bengalensis euptilurus) that lives only on the Tsushima Island in Japan and is threatened with extinction. Because the TLC population is small, genetic management is important. In this study, we obtained the draft genome of the TLC and identified single-nucleotide polymorphism (SNP) markers using a genotyping by random amplicon sequencing-direct (GRAS-Di) analysis. We genotyped 31 captive individuals and 50 wild individuals, of which 48 were from a previous study. The identified SNPs were used to clarify the genetic diversity and genetic structure of the wild and captive populations of the TLC. The size of the genome was estimated to be about 2.42 Gb. The number of SNP markers developed was 139, and although PID and probability of exclusion obtained using these SNP markers were not as high as those reported in the studies of other wild species, these SNP markers could be used to identify individuals and parentage. Moreover, the genetic diversity indices of the captive population were similar to those of the wild population. These SNP markers will be useful for understanding the ecology of the TLC and planning conservation strategies.

The mitochondrial genome of the threatened tideland snail Pirenella pupiformis (Mollusca: Caenogastropoda: Potamididae) determined by shotgun sequencing.

The nearly complete mitochondrial genome of the threatened tideland snail Pirenella pupiformis (Mollusca: Cerithioidea: Potamididae) was determined by shotgun next-generation sequencing. The mitogenome is comprised of 13 protein-coding genes (PCGs), two ribosomal RNA (12S and 16S) genes, and 22 transfer RNA genes (tRNAs). This gene order is consistent with the previously published mitochondrial genomes of other species belonging to the family Potamididae. The family Potamididae including P. pupiformis was recovered as a monophyletic group in the superfamily Cerithioidea.

First mitochondrial genomes of Capitellidae and Opheliidae (Annelida) and their phylogenetic placement.

The complete mitochondrial genomes of Notomastus sp. (15,776 bp) (Annelida: Capitellidae) and Armandia sp. (18,538 bp) (Annelida: Opheliidae) were assembled for the first time. A group II intron (303 bp) was found in cox1 of Notomastus sp. A phylogenetic analysis revealed that Notomastus sp. and Armandia sp. were monophyletic, and this clade was clustered with echiurans, although the possibility of the effect of long-branch attraction should be considered.

Genetic Population Structure of Wild Boars (Sus scrofa) in Fukushima Prefecture.

We aimed to reveal the dispersal and gene flow of the local wild boar (Sus scrofa) population and find their genetic boundary in Fukushima Prefecture. After the nuclear incident in 2011, the land was considered a difficult-to-return zone, and the increase in the number of wild boars was pronounced. To provide an effective management strategy for the wild boar population, we used multiplexed inter-simple sequence repeat genotyping by sequencing (MIG-seq) and clarified the genetic structure of wild boars. We obtained 328 single-nucleotide polymorphisms from 179 samples. STRUCTURE analysis showed that the most likely number of population cluster was K = 2. Molecular analysis of variance showed significant genetic differences between groups of wild boars inhabiting in the east and west across the Abukuma River. The migration rate from the eastern population to the western population is higher than in the reverse case based on BayesAss analysis. Our study indicates that both the Abukuma River and anthropogenic urbanization along the river may affect the migration of wild boars and the population in western was established mainly by the migration from other neighboring prefectures.

Draft genome sequence data of Indian rhinoceros, Rhinoceros unicornis.

The Indian rhinoceros (Rhinoceros unicornis) is a large herbivore found in northern India and southern Nepal. It is a critically endangered species, with an estimated population of approximately 3,600 in the wild. Genetic factors, such as the loss of genetic diversity and the accumulation of deleterious variations, are critical risk factors for the extinction of endangered species, such as the Indian rhinoceros. To support the conservation efforts of the Indian rhinoceros, we assembled its draft genome. The new genomic data will enable the study of functional genes associated with the ecological and physiological characteristics of Indian rhinoceros and help us establish more effective conservation measures. The muscles of an Indian rhinoceros that died from prostration at a zoo were collected, and the samples were stored at the National Institute for Environmental Studies (Tsukuba, Japan). Sequence data were obtained using an Illumina NovaSeq 6000 platform for short reads and an Oxford Nanopore Technologies PromethION for long reads. We generated approximately 235.2 Gbp of data. From these sequences, we assembled a 2,375,051,758 bp genome consisting of 7,615 contigs. The genome data are available from the National Center Biotechnology Information BioProject database under accession number BOSQ00000000.

Potential differences in chitin synthesis ability cause different sensitivities to diflubenzuron among three strains of Daphnia magna.

Ecotoxicity testing of crustaceans using Daphnia magna has been implemented in the chemical management systems of various countries. While the chemical sensitivity of D. magna varies depending on genetically different clonal lineages, the strain used in ecotoxicity tests, including the acute immobilization test (OECD TG202), has not been specified. We hypothesized that comprehensive gene expression profiles could provide useful information on phenotypic differences among strains, including chemical sensitivity. To test this hypothesis, we performed mRNA sequencing on three different strains (NIES, England, and Clone 5) of D. magna under culture conditions. The resulting expression profile of the NIES strain was clearly different compared to the profiles of the other two strains. Gene ontology (GO) enrichment analysis suggested that chitin metabolism was significantly enriched in the NIES strain compared to that in the England strain. Consistent with the GO analysis, evidence of high levels of chitin metabolism in the NIES strain were observed across multiple levels of biological organization, such as expression of chitin synthase genes, chitin content, and chitinase activity, which suggested that the different strains would exhibit different sensitivities to chemicals used to inhibit chitin synthesis. We found that among all strains, the NIES strain was more tolerant to diflubenzuron, a chitin synthesis inhibitor, with a 14-fold difference in the 48 h-EC50 value for the acute immobilization test compared to the England strain. The present study demonstrates that the differences among strains in chitin metabolism may lead to sensitivity difference to diflubenzuron, and serves as a case study of the usefulness of comprehensive gene expression profiles in finding sensitivity differences.